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ABSTRACT. We propose Kolmogorov-Smirnov test as a means for recognising event-by-event fluctuations 
of rapidity distributions in relativistic heavy ion collisions. Such fluctuations may be induced by the spinodal 
decomposition of the rapidly expanding system during the 1st order quark-hadron phase transition. 



1 INTRODUCTION 

In heavy ion collisions matter in state with deconfined 
quarks (quark-gluon plasma) may be produced de- 
pending on the collision energy. The deconfined matter 
expands and cools down very quickly, possibly reach- 
ing the line of the 1st order phase transition between 
the quark phase and the hadronic phase. During this 
transition the system supercools and if the expansion 
rate is fast enough it may go through the so-called spin- 
odal decomposition - the fragmentation of quark-gluon 
plasma into pieces of characteristic size O [TJ [2] . These 
fragments will eventually hadronize, emitting parti- 
cles with a unique rapidity spectrum which reflects a 
unique nature of fragmentation of a given event. Ra- 
pidity spectra thus fluctuate non-statistically on event- 
by-event basis. 

We propose a method for the study of these fluctua- 
tions based on Kolmogorov-Smirnov (KS) test. KS test 
addresses the question if two empirical distributions 
F ni (X), F m (X) are generated from the same underly- 
ing probability distribution. The maximum difference 
D between F ni (X),F n2 (X) is distributed according to 
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with n = n\n<2,j[n\ + n-z). 

2 RESULTS 

Below we calculate quantity Q = 1 — H for rapid- 
ity spectra (empirical distributions F ni {X), F n2 (X)) of 
random pairs of simulated events where rapidities were 
first generated according to test (Gaussian) distribu- 
tions and then for pairs of events generated with our 
Monte Carlo droplet generator which simulates particle 
emission from the fragments emerging from spinodal 
decomposition. If rapidity spectrum for each event is 
generated from the same underlying distribution, we 
expect an equal number of event pairs at any Q - a 
uniform, constant line histogram (we expect this in the 
case when no fragments are present). For events drawn 
from different underlying distributions the pairs group 



This point is illustrated in Fig. 1 where we show 
results of KS test on data randomly generated from 
two different Gaussian distributions: a) half of 10 5 
events was generated according to the distribution 
(fj.,a) = (1,0.1), half according to (/x, a) = (2,0.1), 
b) half according to (/i, a) — (1.99,0.1), half accord- 
ing to (fJ.,(j) = (2,0.1). Multiplicities ranged in both 
cases from 8 to 535. The Q value was calculated for 
10 5 random pairs of these events and the histogram of 
the number of pairs as a function of the corresponding 
Q value was plotted. 



KS prob, V=20 lm"3. Gauss! = (1,0.1), Gauss2 = (2,0.1) U 10"5 random pairs/10' 
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K5 prob, V=20fm"3, Gaussl = (1.99,0.1), Gauss2 = (2,0.1) !, 10*5 random pairs/1 0*5 events 
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Fig. 1. KS test for data from two different Gaussian 
distributions: a) (/i, cr) = (1,0.1) and (/x, a) = (2,0.1), 
b) (n,a) = (1.99,0.1) and (/*,cr) = (2,0.1). 

We can see in Fig. la that half of pairs is distributed 
uniformly (5000 pairs in each bin) and half (50 000) is 
concentrated in the first bin. The first half consists 
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events) which form a constant line, the second half 
consists of events from opposite distributions (unlike 
events) which group at small Q values. Note that the 
statistical fluctuation is ±y/N. The deviation from the 
constant line at N = 10000 above y/N = 100 repre- 
sents non-statistical fluctuation. 

Fig. lb shows the power of the KS test. Although 
the two Gauss distributions differ by just 0.01 in their 
mean value (a = 0.1), the KS test clearly recognized 
that the events were not drawn from a single distribu- 
tion. 

After the test of KS test on two Gaussian distri- 
butions, we applied the KS test on data generated 
with the Monte Carlo droplet generator. This gener- 
ator simulates particle rapidity spectra emerging from 
the fragments (droplets) originating in the spinodal de- 
composition. Droplets are generated from a blast-wave 
source with tunable parameters. In Fig. 2 we show re- 
sults of KS test on droplets of size a) 50 fm 3 and b) 10 
fm 3 at temperature T = 175 MeV. As before, we used 
a sample of 10 5 event pairs randomly chosen out of 
10 5 events. All particles come from the droplets. Non- 
statistical fluctuations are clearly visible in both cases 
and grow with the droplet size. 



KS prob. V=10fm»3, 20% from droplets, T= 175 MeV, 10 n S random pairs 10*5 events 
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V=10 fm"3, T = 175 MeV, 10 ft 5 random pairs/10 A 5 events "3 
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Fig. 2. KS test on droplets of size a) 50 fm 3 and b) 
10 fm 3 at temperature T = 175 MeV. 

Finally, we fix the droplet size at 10 fm 3 (T = 175 
MeV) and vary the percentage of particles originating 
in the droplets. In Fig. 3a 20% of particles comes from 
the droplets as compared with 60% in Fig. 3b. Again 
the sample was 10 5 event pairs chosen randomly from 
10 5 simulated events. The non-statistical fluctuations 
are present even for the 20% case. They vanish as ex- 



24000 
22000 
20000 
18000 
16000 
14000 
12000 
10000 
8000 



h3 

Entries 99998 
Mean 0.3995 
RMS 0.31 



(a) 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

KS prob, V=10 fm"3, 60% from droplets, T= 175 MeV, 10"5 random pairs/10"5 events | 
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Fig. 3. KS test for the case when only a fraction of 
particles comes from the droplets: a) 20% and b) 60%. 



3 CONCLUSIONS 

We studied event-by-event fluctuations of rapidity dis- 
tributions by means of the KS test. We have shown 
that the KS test can recognize whether or not the 
rapidity spectra are drawn from the same underlying 
distribution. If spinodal decomposition occurs during 
the 1st order quark-hadron phase transition, the result- 
ing droplets lead to many underlying distributions and 
thus to event-by-event fluctuations. As a next step, we 
plan to apply this method on data by NA49 collabora- 
tion. Of course, if the fluctuations are indeed there, it 
does not follow automatically that they are caused by 
the droplets. 

This research has been supported by VEGA 
1/4012/07. 

REFERENCES 

[1] B. Tomasik et al, contribution to this conference, 
arXiv: 0711.4932 [nucl-th]. 

[2] O. Scavenius et al, Phys. Rev. D 63, 116003 
(1993). 

[3] I. N. Mishustin, Phys. Rev. Lett. 82, 4779 (1999). 



